Can an evolving Universe host a static event horizon? 
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We prove the existence of general relativistic perfect fluid black hole solutions, and demonstrate 
the phenomenon for the P = wp class of equations of state. While admitting a local time-like 
Killing vector on the event horizon itself, the various black hole configurations are necessarily time 
dependent (thereby avoiding a well known no-go theorem) away from the horizon. Consistently, 
Hawking's imaginary time periodicity is globally manifest on the entire spacetime manifold. 



INTRODUCTION 



There is by now abundant evidence that our Universe 
is expanding on average. Nevertheless, much of the vi^orks 
on black holes, which play a central role in modern astro- 
physics, focus on stationary and asymptotically flat situ- 
ation. It is therefore desirable to have black hole models 
embedded in cosmological environment to see if those 
objects can reveal some unexpected features of the un- 
derlying theory of gravity. Historically, an investigation 
of the effects of the cosmological expansion on local sys- 
tems, was motivated by the question of whether an atom, 
a star. Solar System, galaxy or any other bounded sys- 
tem expand following the rest of the universe. Although 
this question has a long history dating back to the 1933 
paper by McVittie [l] introducing a spacetime metric 
that represents a point mass embedded in a Friedmann- 
Lemaitre-Robertson- Walker (FLRW) universe, one still 
lacks an affirmative answer to this open problem in gen- 
eral relativity. In fact, the physical properties of the 
McVittie solution are an active field of research and dis- 
cussed by several authors, e.g^ 0, Q, [Ij, 0. Later 
work by Einstein and Straus \i\ describes a patchwork of 
Schwarzschild black holes with an FLRW universe, while 
subsequent generalizations of their ideas [1| replace the 
external FLRW metric by the inhomogeneous Lemaitre- 
Tolman-Bondi (LTB) dust time-dependent solutioii [a 
and by its non spherically-symmetric generalization [lOj . 
The former, however, are time symmetric black hole so- 
lutions and so they do not describe at a satisfactory level 
any phenomenon associated with a dynamical black hole 
in an evolving Universe. Such processes are expected to 
play a significant role in black hole formation, gravita- 
tional collapse, evolution of primordial black holes, etc. 
This has invoked numerous studies, e.g., 11 1, 12 1, [l3| . 
[l3 |. [isl where new non-static solutions have been con- 
sidered. Overall the references to the inhomogeneous 
solutions in general relativity are quite extensive, e.g., 
see [l3|, [l3] and references within. In this work we are 
mainly motivated by the question if an evolving and non- 
stationary universe can host a static event horizon. 



The event horizon is a central object in the descrip- 
tion of a black hole, with the Schwarzschild solution be- 
ing its most prominent representative. The near hori- 
zon structure of the latter is a product of a Rindler 
space and a two sphere, which we write as Rindler x S^, 
and defines the corresponding causal structure. Specif- 
ically, the Rindler piece in the metric typically char- 
acterizes the event horizon as a boundary of region in 
spacetime from behind which no causal signals can reach 
the observers sitting far away at infinity. The Euclidean 
Rindler metric that is obtained by substituting t = —He 
yields fiat two-dimensional Euclidean metric written in 
polar coordinates, provided the angular variable has the 
correct periodicity. If the periodicity is different, then 
the geometry would have a conical singularity at the 
would have been horizon x = 0. Finally, the relation 
/? = HAt between the periodicity At in Euclidean time 
tE, and the inverse temperature /3, leads towards the 
celebrated Hawking-Bekenstein black hole temperature. 
Even though we discussed the temperature in the context 
of vacuum Schwarzschild black hole, the concept of black 
hole temperature that owes its life to the "no conical 
singularity" demand, is true even in a presence of matter 
because the near horizon geometry is always Rindler-like. 
In this work we take one step further and assume that a 
horizon and its underlying Rindler-like structure in the 
near-horizon limit, also exists in the time dependent and 
spherically symmetric case. In particular we will assume 
that on the horizon itself the metric is static, i.e. we 
postulate the existence of a local time-like Killing field. 
This in contrast to the Schwarzschild, or Schwarzschild- 
deSitter black holes where time-like Killing fields are 
manifest in any point. 

The rest of this work is organized as follows. First, we 
shortly discuss the McVittie solution, analyze the gen- 
eral structure of the perfect fiuid energy-momentum ten- 
sor, the emerging baryon number conservation and the 
Euler relativistic equations. Then we introduce the per- 
fect fluid isotropic condition and the perfect fiuid equa- 
tion of state which constitute our basic equations for the 
corresponding spherically symmetric metric components. 
Their solution yields the corresponding spherically sym- 
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metric configuration in the presence of a perfect fluid. 
We solve those equations perturbatively, order by order, 
with dimensionless radial coordinate serving as an expan- 
sion parameter, around the local Killing horizon. The 
Rindler metric on this static surface serves as boundary 
condition in our iterative series solution, and defines the 
corresponding causal structure. 

The McVittie solution 

One of the most prominent representatives of non- 
homogeneous and non-stationary exact solutions of Ein- 
stein field equations (adopting from this point onward 
the natural units 16ttGn = c = 1) 

is the McVittie line element, in its special fc = case, 
explicitly given by 

Here, /i is a dimensionless quantity defined by 

r is the radial coordinate and dfl^ = dO^ + sm{9)'^ dcj)^ is 
the solid angle on a 2-sphere. The corresponding metric 
is an exact solution of Einstein field equations for any 
constant m, provided the corresponding energy density 
is a function of the time coordinate alone. It approaches 
the FLRW universe in the /i <C 1 limit, and tends to 
the Schwarzschild solution in the static case. Although 
admitting those two tenable limit cases, it does not cover 
other cases of interest in the interface between cosmology 
and local black hole physics, and cannot be regarded as 
general problem solution for the following reasons: 

(i) Although McVittie's solution tends to the pure 
FLRW solution in the far region, it falls short to match 
the standard description in terms of a perfect fluid in 
the near region. Specifically, any barotropic equation of 
state P = P{p) where p and P are the energy density and 
isotropic pressure as measured in the instantaneous rest 
frame of the fluid at the corresponding point, is incom- 
patible with this solution. In particular, assuming that 
the underlying matter is described by the perfect fluid 
energy-momentum tensor given by 

T^"" ^-{p + P)U''U'' -Pgf"" (4) 

where If^ is the four velocity vector that satisfies the 
normalization condition 



one is led towards p which is a function of time <, implying 
that there is no -P = P(p) equation of state. 

(ii) McVittie's line element is over restrictive since it 
prevents the hole from accreting, and enforces its energy 
density to be exactly homogeneous on some set of spatial 
slices. 

Additional analysis of McVittie solution may be per- 
formed by using the concept of the apparent horizon 18]. 
The latter defined as a locus of vanishing geodesic ex- 
pansion, which in the context of the McVittie's metric 
H yields 

1 _ ^ _ = 0, (6) 

?' 

which expresses r as a function of t through the time- 
dependence of the Hubble parameter H. As time goes 
on the area of the apparent horizon changes, in such a 
way so that r„-the smaller root of Eq.® moves inward 
as H decreases with time towards its asymptotic de Sitter 
value, while the larger root moves outward. In the limit 
t — )■ oo as the Hubble parameter tends to some asymp- 
totic positive value H — > Hq > 0, one may show that 
r = r_, < = (X) is a regular event horizon. In the partic- 
ular Hq = case, on the other hand, the r — r_ surface 
does not acquire a straightforward meaning as an event 
horizon, due to soft singularities of scalars constructed 
from quantities that involve two derivatives of the Rie- 
mann tensor. In this work we deal with a static event 
horizon. 

To summarize, although the McVittie solution de- 
scribes a black hole, at least for positive Ho, it still does 
not encompass some effects and properties one would ex- 
pect for a physical black hole in a universe full of matter 
or radiation. 



Setup and Field Equations 

We now describe our basic assumptions and main equa- 
tions that we follow below. In this work we consider 
a time-dependent and spherically-symmetric solution of 
Einstein field equations Eq.(IT]) in the presence of per- 
fect fluid Eq. dU . We adopt the isotropic coordinates and 
consider, without loss of generality, the following line el- 
ement 

ds^ = -T{t, r)de + R{t, r) {dr^ + r^dil^) . (7) 

where T(t, r) and R{t, r) are functions of the time t and 
the radial coordinate r. Introducing a dimensionless vari- 
able X and the characteristic mass scale m via 



U^U^ = -1, 



(5) 



r = -{l+x) 



(8) 
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we may rewrite the line element Eq.© as 

w?R{t, X 



-T{t, x)dr 



l + x) dVi' 



(9) 

where for simplicity we still follow the same notation for 
the metric components as in Eq. ([7]) . Recognizing the fact 
that near Schwarzschild black hole event horizon the met- 
ric acquires a typical Rindler x structure, we assume 
that as we approach towards the surface a; = the metric 
tends to 



ds^ 



-dt^ + An? {dx^ + d^"^) 



(10) 



thereby capturing the causal structure associated with 
the corresponding Rindler line element in the close vicin- 
ity to the a; = event horizon. We can easily convince 
ourselves that the constant m coincides with the mass in 
the Schwarzschild solution. The latter is obtained by re- 
placing a{t) in Eq.® with unity and tends to Eq. pUl) in 
the a; — >■ limit. As we see below, in a more general cases 
the constant m does not acquire such simple meaning. 

The line element Eg. dTUl) has no conical defects in its 
Euclidean regime, once we have assumed that its Eu- 
clidean time coordinate has the correct periodicity At 
given by 



At = 87rm. 



(11) 



Formally, the metric given by Eq. ((TO)) serves as a bound- 
ary condition and also as the zeroth order in our series 
solution 

2 

nt, - ^(1 + xh{t) + x^hit) + . . . ) ^^2) 

Rit, x) = 16 (1-1- xgi(t) + x^g2{t) + ...), 

where /„(i),(;„(t) (n — 1,2,...) are yet undetermined 
functions (to be referred below as harmonics). An ex- 
pansion around black hole horizon, in a non-cosmological 
context, has been performed in 23[. Interestingly, as we 
show below, under the assumption that the line element 
Eq. (jlOp on the horizon has no conic singularity, i.e. satis- 
fies Eq. (Ilip , the subsequent terms in the expansion must 
respect no conic singularity as well. In other words, as- 
suming Eq. (jll[) for n = terms is sufficient to make all 
other n > terms periodic in the Euclidean time coordi- 
nate. 

We now derive our main equations given by Eq. (jl7p 
and Eg. dT^ below, utilized to obtain the analytic ex- 
pressions for the functions fn{t), gn{t) which were already 
defined in the metric expansion Eg. dT^ . The obtained 
solution corresponds to a time-dependent and spherically 
symmetric radial flow of a perfect fluid, such that the 
four velocity vector has a non- vanishing radial compo- 
nent U"^ ^ 0, and may be parametrized according to 



cosh((^(t, r)) smh{Lp(t,r)) 



0,0 



(13) 



once the rapidity ip{t^ r) has been introduced. Plugging 
the latter into Eq. (j4|) one obtains an explicit expressions 
for the corresponding components of the perfect fluid 
energy-momentum tensor, expressed through the follow- 
ing useful combinations A, B, C, -D, defined by 

A = -2gl ^Tt ^ p cosh^ {(p{t, r)) + Psinh^ {(fit, r)) 
B = 2g;. = -T; = psmh\^{t,r))+Pcosh'{^it,r)) 

c - 2g« = 2gt - -t£ = -rf - p 



D = -2, 



R{t,ry \l R{t,r) 
{p + P) smh{(f{t^ r)) cosh((y9(t, r)). 



(14) 



Here, we also utilized Einstein equations, i.e., propor- 
tionality between the Einstein tensor and the energy- 
momentum tensor ■ The specific algebraic structure of 
the perfect fluid energy momentum tensor, implies that 
the emerging combinations A, B, C, D are not indepen- 
dent. In fact, they are related by 



AB-D^ = C{A -B + C). 



(15) 



Now, assuming the spacetime metric admits the form 
given by Eq.®, one views Eq. (fT5|) . as an equation for the 
metric components T(t,x) and R{t,x). Formally speak- 
ing, by eliminating p, P, and ip we have inverted the equa- 
tions Eq. (jl4p and then by using the Einstein equations we 
ended up with the so-called perfect fluid isotropic condi- 
tion. This condition was originally obtained by McVittie 
and Wiltshire [l^ , and expresses the necessary condition 
on the metric to describe spherically symmetric flow of a 
perfect fluid. Other works where time dependent spheri- 
cally symmetric configuration have been studied may be 
found at We should mention that in case 

we insist on the comoving solutions, which in principle 
are always possible for a single-component perfect fluid, 
Eq.dlSl) is solved for = and = which lead 
to J7'' = 0. In this work, however, we consider solutions 
which are easier to obtain and study in non-comoving 
coordinates. Numerous works, including the LTB dust 
solution 0, where the comoving coordinates are consid- 
ered may be found in the literature. 

The perfect fluid isotropic condition is supplemented 
by matter equation of state, which in this work is chosen 
as the linear barotropic 



P = wp 



(16) 



where w is some constant usually bounded between plus 
and minus unity. By utilizing the combinations that were 
defined in Eq. , the equation of state Eq. ([T5)) may be 
rewritten as 



C = wiA- B + C), 



(17) 
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while the isotropic condition Ea. (jl5|) acquires a simpler 
form 



^2 

AB-D-' = — , 
w 



(18) 



for It; 7^ 0. Clearly, the w — case should be tackled with 
Eg. (ITS]) . Together the equations (|17l) . (IT51) constitute a 
system of two basic coupled equations for the two un- 
known functions T{t, x) and R (i, x) which describe radial 
perfect fluid flow subject to the corresponding equation 
of state. Those equations are difficult to solve even nu- 
merically and in the remainder of the paper we will solve 
them perturbatively order by order around the local hori- 
zon, according to the scheme mentioned near Ea. (ll2p . 

Finally, let us consider the underlying dynamics of 
the perfect fluid matter, which is governed by the in- 
herent Bianchi identity and the corresponding energy- 
momentum conservation equation 



V^T'^'^ = 0. 



(19) 



In case of perfect fluid energy-momentum tensor Eq.Q 
the latter may be re-written as follows 



(p + P) (C/'^V^C/^ + U^^Vt^U") 
{p + P) WW + Vt^Pg^" = 0. 



(20) 



The particular form of the energy-momentum tensor 
Eq.Q, implies that its divergence, given explicitly by 
Eq. ipU)) is a sum of two terms. Specifically, one term is 
parallel to the four velocity vector while the other is nor- 
mal to it. In fact, projecting Eq. ([20l) on Ui, we find the 
component that is parallel to the vector 



[pW) + PV^W - 0, 



(21) 



which is known as the baryon number conservation equa- 
tion. The other term, known as the relativistic Euler 
equation, is orthogonal to the four- velocity and is explic- 
itly given by 

(p -I- P) WW^U" + V^P {WW + cf) = 0. (22) 

Let us notice that Eq. ([2T|) can be written as a vanishing 
divergence of the current 



(23) 



once we have introduced the proper density n (number 
of particles per volume) 



3' 



according to 



nW 



r 



(24) 



(25) 



In the particular case of a linear barotropic equation of 
state given by Eq. (|16p . the proper density n acquires the 
following simpler form 

n = p^. (26) 



As a final remark we note that once the metric and the 
corresponding Einstein tensor have been worked out, the 
following expressions for the rapidity ^p{t,r) 



sinh^((^(i, r)) 
cosh^ {if {t,r)) - 



B-C 



A- B + 2C 
A + C 



(27) 



A-B + 2C 

may be utilized to determine the perfect fluid radial flow. 

NO-GO THEOREM FOR THE STATIC CASE 

We now prove that the only perfect fluid P = wp 
static, spherically symmetric black hole solution is the 
Schwarzschild solution with vanishing P and p. For sim- 
plicity, we discuss flrst the w — case and turn to w ^ 
case afterwards. Then, we also consider the static case 
limit which proves to be useful when compared to the 
time-dependent case in the following. 

The lii = case 

Let us flrst consider the simpler static case with a van- 
ishing proper pressure P = 0, which corresponds to the 
w = case. Assuming the spherically symmetric line el- 
ement given by Eq.Q let us deflne for convenience the 
functions f{r),g{r) 

T{r) = e-^^'') , R{r) = e^^'') (28) 

and their derivatives with respect to the radial marker r 

fir) = F{r) , g'{r) = G(r). (29) 

Substituting now the corresponding metric components 
into our basic equations Eq. lfTSl) . Eg. dTTl) we are driven 
towards the following flrst order differential equations for 
F{r) and G{r) 



2F{2 + rF) + 2{G + r{F' + G')) = 
2P(2 + rG) + G{A + rG) = 0. 

Solving those gives rise to 

Srn 



(30) 



F{r) 
G{r) 



4to 



(31) 



r(m + 2r) 

and the corresponding Schwarzschild line element 



"■°Sch — 




m \ 



1 + 77^) (dr^+r^dn^) 
2r J 



(32) 
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This proves that the only static, spherically symmet- 
ric solution with dust perfect fluid is the dustless 
Schwarzschild solution. 

Notice that in this simple case the proof did not rely 
on the assumption that local Killing horizon exists. In 
the following w ^ Q case, however, we are not able to 
prove our theorem without assuming the existence of a 
local Killing horizon. 



From this expansion we can learn that it is w indepen- 
dent and thus describes a Schwarzschild black hole. To 
summarize, we have proven that any static metric with a 
Killing horizon in the presence of a perfect fluid, is nec- 
essarily a Schwarzschild solution with a vanishing proper 
energy density and a vanishing proper pressure. Inter- 
estingly, in the w = case, considered before, this claim 
holds even if we do not assume the Killing horizon exists 
at X = 0. 



The w 7^ case 

Now let us consider the more interesting w ^ static 
case, with a non- vanishing pressure. In such case our 
basic equations Eq.(fT7|). Eq. (fT8|) are more involved, and 
we are not able to solve our basic equations analytically. 
Nevertheless, assuming the existence of a local Killing 
horizon aX x — dictates Rindler x S^ metric given by 
Eq. pO)) as a; — > 0. This fixes the zcroth order terms in the 
expansion Eq. (jl2p and proves to be useful for a construc- 
tion of our solution as a power series in x, around x — Q. 
Specifically, we assume that the functions F{r),G{r) de- 
fined above, admit the following expansion 

F — \- an + aix + a2X'' + . . . 

X (33) 

G = Po + l3ix + hx^ + . . . 

where a„_i and /?„ (n = 0, 1, 2, . . . ) are the correspond- 
ing constants that we find order by order by plugging the 
series Eq. ipS)) into our basic equations Eq. lfTT]) . Eq. (fTSl) . 
The latter may be written as 



(34) 



P(i)(r)+ii;P(2)(r) = 
p2)(r)-Hu;P(3)(r)-0, 

respectively, where (r), P(2) ^(3) (^) are w inde- 
pendent quantities that are defined according to 

P(i) (r) = (2 + rF)F + 2(G + r(P' + G')) 

P(2) (r) = (2 - rF)F + lOG + 2r(PG + G^ - F') + G' 

P(3)(r) ^ {A{F + G) + 2rFG + rG^){m + rG^ +ArG'), 

(35) 

where the prime denotes derivative with respect to r. We 
can straightforwardly show that the previous solution in 
the ui = case, is still a solution of Eq.(|34|) for any 
ui 7^ 0. Furthermore, using the small x expansion we 
can verify that in fact this is the only solution, implying 
that P(i) (r), P(2) (r) and P(3) (r) all vanish. Having solved 
Ea. (l34l) for the lowest order, i.e. a_i and /3o, we proceed 
to the next, and so on solving each time for the leading 
terms. Following this procedure we are able to solve for 
all the following orders, and we bring here some of the 
terms 



1 



a-i 



-;ao 



m 



-;/3i 



1 

'2m 
3 

2m 



1 



4m 
7 

4m 



;Q!2 



2 

m 



(36) 



Perfect Fluid in a Static Black Hole Background 

Now let us consider the perfect fluid flow in the static 
case limit 2J|, which provides some useful insights into 
the time-dependent case that follows afterwards. In the 
static case the relativistic Euler equation Eq.([22|) is char- 
acterized by the time independent metric components 
T{t, x) = T{x) and R{t, x) = R{x) and by the fact that 
the fluid is static. This dictates that the four-velocity 
vector C/'^ is purely timelike 



f/° = i-gooY 



= 



for A 7^ (37) 



and also leads to the following identity 



0"" 2^ dx-- 



(38) 



Utilizing now Eg. (1571) and Eq. ([55|) one may rewrite the 
relativistic Euler equation given by Eq. ((22|) . for the space 
components v — i {i — 1,2,3), as following 



P. 



T 



-T, 



(39) 



(for V — Q the relativistic Euler equation is satisfied triv- 
ially in the static case). The latter is nothing but the 
non-relativistic equation for hydrostatic equilibrium in 
static gravitational field and has few interesting impli- 
cations. In particular, assuming that such a spacetime 
is also inhabited by a static black hole, then its event 
horizon (which also coincides with its Killing horizon in 
this static case) is identified by Tj^^^ =0 where ru stands 
for the horizon's coordinate. We notice now that under 
the assumption that (p + P)|r,, 7^ and finite Tilr^ on 
the horizon rh. Eg. (1391) implies that the pressure gradi- 
ent is necessarily divergent. Therefore, avoiding this un- 
physical divergence calls for a vanishing of the following 
combination 



(40) 



on the horizon rh- In case we also insist on the equation of 
state Eq. (fTB)) . then we are also driven towards a vanishing 
proper energy density on the horizon 



pV 



0. 



(41) 
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As we will see below, the severe restriction given by 
Ea. (l40l) is elegantly removed once the time dependence 
is introduced. 



only p 7^ 0, g 7^ can go beyond. In fact, calculating 
the proper energy density p on the horizon up to the first 
order, yields 



DYNAMICS: HARMONIC EXPANSION 



P = 



1 + w 



0{x) 



(48) 



First Harmonic 

We now turn to discuss the time-dependent solutions. 
As we have mentioned above we hereby assume that our 
spacetime hosts a local Killing horizon, implying that as 
X — >■ the components of our metric Eq.([9]) admit the 
expansion given by Eq. ()12p . In particular, in the leading 
order the corresponding line element is given by Eq. (flUl) . 
The governing equations for the first harmonics /i {t) and 
(7i(t), are found once we plug the expansion Eq. (jl2l) into 
our basic equations Eq. (fT7|) . Eq.lfTS]) and keep the leading 
order in each one accordingly. This way, we obtain the 
following equations for the first harmonics 

1 + w + (1 - w) /i (t) + wgi [t) - mm^g'l {t) ^ (42) 



that depends on the product pq. This just refiects the 
fact, that our solution describes a static horizon, and 
therefore any shift of the time coordinate t t + At 



results in p — > pe" , q qe^'^ so that the product 
pq — > pq stays unchanged. The non-vanishing proper en- 
ergy density on the horizon, in this time-dependent case 
is in contrast with our previous comment near Ea. (j41l) . 
where we have indicated that the proper density p should 
vanish in a static case. 

Interestingly, the first harmonic may be fixed by an 
alternative demand that p and P are non-singular on 
the horizon, without providing any particular equation 
of state. To this end, let us calculate to the leading order 
the proper energy density, the pressure, and the perfect 
fluid isotropic condition Eq. p^ 



1 + /i (t) - IQm'g'l it) = 



(43) 



which are valid for any choice of the parameter w. The 
latter are linear and admit an analytical solution explic- 
itly given by 



5i(0 = -l+pe"* + ge-"* 
/i (t) 



2 + pe'^* + ge-"* 



(44) 



where the Euclidean angular frequency w has been intro- 
duced according to 



l/4m 



(45) 



Here we notice, that in order that our expansion stays 
valid it is actually needed that the following product re- 
mains small 



(46) 



From the first order terms Eq. (j44p given above, we can 
already learn about one of the major features of our time- 
dependent solution which holds to higher harmonics as 
well. Particularly, it has no conical singularity once con- 
tinued to the Euclidean time via t — —itE- In fact, keep- 
ing in mind that the correct periodicity Eq. (ITlT) of the Eu- 
clidean metric Ea. dTO]) is assumed to hold, we can check 
by substitution that a shift At, 



± iuj{tE + At) = ±iujtE ± 27ri 



(47) 



gives rise to the 27ri shift. 

Dropping the time-dependence by considering p = q = 
case simply brings us back to the Schwarzschild, and 



P = 



Au' {~1 + hit) ~ g,it)) 



4c.2(l + A(t)-^<(i)) 



-0{x") 



(49) 



+ 0{x°) (50) 



AB-D"^- C{A -B + C) 

(-1 + f,{t) - g,it)) (1 + A (t) - ^g'{{t)) 



+ 0(.T-i), 



(51) 



respectively. Assuming now the leading terms in Eq. (j49p 
and Eq. ([5ni) vanish, we are guaranteed that perfect fluid 
isotropic condition is satisfied. The perfect fluid is non 
static by construction, see Eq. (fT3l) . Furthermore, Eq.([27|) 
then tells us that associated with our solution is the ra- 
pidity 



sinh(2(y5 {t, r)) 



2pq 



+ Oix). (52) 



The important point is not that T^i/ is non-diagonal, but 
that there is a radial current or flow represented by this 
rapidity. 

But this is not all there is; the physical quantities such 
as proper energy density p are related to the Ricci tensor, 
which involves second order metric derivatives. There- 
fore, it seems that we must go even further, at least to 
the second harmonic for additional contribution. 
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Second harmonic 

Before we consider the full set of equations for the 
second order terms f2{t) and 52(^)1 let us try to learn 
how the restriction Eq.pO]) is removed once the time- 
dependence is introduced. To this end, we consider the 
expansion Eq. p^ up to the second order, with the first 
harmonic Ea. (l44p plugged in. This way, we can learn 
about the role of just the perfect fluid constraint Eg . (fTSl) . 
but still without specifying an equation of state. Of par- 
ticular interest are the proper energy density, and the 
proper pressure that may be expressed up to the second 
order terms, according to 



pe^' -qe-^' -pq + f2{t)-92{t)] (53) 



we imposed also P = wp we would have faced another 
full differential equation. In fact, the relevant equations 
for 52(0 and f2if) are obtained by substituting the near 
horizon expansion Eq. ()12|) with already known first order 
terms Eq.(|ii|). into our basic equations Eq. pT)) . Ea.([T^. 
Keeping the leading terms up to the second order we de- 
rive the equations for the corresponding /2(t) and g2{i)- 
Quite generally, Eq. ([T7|) does not comprise any deriva- 
tives of the function /2(i), which allows to express the 
latter through 52 (i) and its derivatives according to 



f2{t) = 



1 



12(-l + w) 
6g(3 + 2w)e""* 



- 23 - 33a; + 6p(3 + 2w)e'^ 



+ 4(1 + 37«)g2W- 4.92 W) 



I2pq{-1 + w) 
4 



(58) 



P = Auj 



21 Z^ + ^^^t 



9g 



ipq + Sf2it) + g2{t)- —g'^it) 



(54) 



Keeping in mind that the sum p + P is significant in the 
static case, we bring here the corresponding expression 
up to the second order 



p + P = 2uj^(5 + Spe'^* + 3ge-'^* - 12^^ 



l2f2{t)-4g2{t)-—gUt) 



(55) 



and also the static case limit of the expressions given by 
Eq.®, Eq.dSl, Eq.dSS]) 



p ^ 3^^n + 4f2 - ig2) 
P = w2 (-23 +12/2 + 4.92) 
p + P = 2ij^{5 + 12/2 - 432) 



(56) 



Interestingly, there is a close relation between the perfect 
fluid isotropic condition Eq. (IT51) and the vanishing of the 
mentioned combination p + P = 0. Specifically, as may 
be verified with the help of MathTensor package (some of 
the expressions are too long to be written here explicitly) 
the perfect fluid isotropic condition is equivalent to the 
condition (p + P)\r,^ = in the static case. In such case 
characterized by /2(^) = /2 and 52 (^) — 92, where /2 and 
g2 are some constants, those two conditions imply 



92 



+ 3/2 , p=-P = 6LoH3-Af2), 



(57) 



still without specifying a concrete equation of state. An 
attempt to impose on top of it an equation of state of 
the type P = wp, will take us back to Ea.p6| and to 



the corresponding coefficients given /2 = f and (?2 = |, 



which brings us back to Eg. (j4T|) . 

Once t— dependence is introduced, imposing p + P — 
leaves us with a full differential equation for 32 (0- Had 



Plugging this relation into the corresponding leading or- 
der of the perfect fluid isotropic condition Ea. ([T^ . we are 
led towards the following differential equation for g2{t) 



-ws"{tf + 2a;2(l + wfs{t)s"{t)- 
t^2(l - wfs'{tf - Auj\l + wfs{tf = 



(59) 



where 



s{t) =14 - 10 (pe"* -f qe"'^*) 

*)-452(t). 



3 (p2e2"* + q2e-2-^^ 



(60) 



In a particular case of a dust perfect fluid, associated 
with It; = 0, the latter can be furthermore simplified 



The solution of Eq. ([5^ is evidently given by 



(61) 



s{t) = -4ao (2(1 -w) + {\ + w)2e2-(*-/5o) + 5-2^(4-/30)^ 

(62) 

which according to Eq.lfSS]). Eq. (15(11) yields an explicit 
expressions for the second order terms 



,. / N 3 ^ 1 + 3w 
Ut) = - + 1 - £— 



2 2 

P 2uit , q -2ujt 

—e + e 

77 1 ~ ^7 



92{t) 



i,pq 



3(1 - w) 5 
2{l + w) ^ 2 



qe 



qe 



2 2 

p g2wt _j_ q Q-2'^t 

1-7] 



(63) 



where the constants k, "q, and ^ have been introduced 
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according to 



2ao = -^^pq 



(64) 



4ao 



1-rj 



1 



The coefficients k^r],^ are not independent and in fact, 
are found to respect the following consistency relation 



^2^^ (l-fc)fc 



(1 - 77)77 ■ 



(65) 



We now discuss the physical quantities, such as the 
energy density and the rapidity associated with our so- 
lution. In fact, the latter is non-trivial since we allow an 
accreting solution with non-vanishing radial component 
of the four- velocity vector, as specified by Eq. (fT3| . The 
proper energy density on the horizon and the rapidity are 
given now by 



P = 



l + w 



0{x) 



sinh(2v?(t,r)) = 

q^rj (1 - fc - 77) e"^'^* + p'^ {r] ~ k) (rj - 1) 
2^pqi] (77-1) 



(66) 



Oix), 
(67) 



respectively. As expected, since the proper density and 
the rapidity depend on the second derivatives of the met- 
ric, the obtained relations Eq. (pS|) and Eq. (p7|) differ from 
the previously obtained expressions Eq. pS)) and Eg. dS^ . 

From expressions Ea. ([M)) and Eg. ([57)1 we learn that: 
(i) The local Killing horizon is non-singular at cc = for 
any constant w (except w ~ —1). (ii) The local Killing 
horizon hosts non vanishing energy density, (iii) The 
combination which appears to be of physical significance 
is ^pq and not just pq. While (i) and (ii) are evident 
from Eq. (|B5| . (iii) has additional implications and de- 
serves further discussion. In fact, as we demonstrate be- 
low, various scalar invariants are all proportional to ^pg 
implying also that the physical quantities p and P are 
proportional to this combination as well. This suggests 
that the expression pe"' -I- ge~'^* which appears in Eq. p^ 
can in fact be gauged away completely, i.e. p,q 0, 
provided we insist ^pq — > const. Specifically, after taking 
this limit, the first order term given by Eq. (|44|) turns into 



/i (t) = -1 
91 W = -2. 



(68) 



In the second harmonic, given by Eq. (p5|) . taking the 
product p, g — > but insisting on keeping the product 



£_pq constant we are driven towards ^ ~ — , that should 
agree with the consistency relation Eq. (|S5]) between ^, p 
and q. Under the following redefinition 



pq' 



1-77 



77 y p"' 

fc2 



1 



77(1 — 77) ri {\ — -q) p'^q^ 

the redefined consistency relation in the 7), q 
reads 

^2 



(1 - r/) 77 



(69) 



limit 



(70) 



for < 77 < 1. Furthermore, utilizing the freedom of 
choice of coordinates 1-^1 + Ai, that does not change ^, 
without loss of generality we can always choose ie^""^' = 

This way the second harmonic may acquire 



_fe_„-2wAt 



the following form 



92 it) 



3(1 



2 (l + w) 
while the rapidity is given by 



-1^ cosh {2ujt) 
3 

cosh {2ujt) 



(71) 



sinh(2(p(t, r)) = - sinh {2ujt) + 0{x^). 



Third Harmonic and Curvature Scalars 



(72) 



We may follow the general prescription to the third 
{n = 3) order. Specifically, let us substitute the metric 
components given by Eq. (|12p , such that the only present 
terms are of order three and lower. Keeping in mind that 
the lower orders n = 1 and n — 2 are already known and 
are given by Eq. (|68| and Eq. ((7T1) . respectively, we expect 
to obtain equations for fsit) and 53 (i). Similarly to what 
we have done in the lower orders, Eq.lfTTl) allows to ex- 
press 53 (t) through fzit) and its derivatives. Plugging 
the resulting expression for 53 (t) into our second basic 
Eq. (|18p we end up with a second order differential equa- 
tion for fsit). Solving this equation for /^{t) yields 



f3{t) 



1 l + 3w 

2 ^ 1 + 77; 



i(l + 37«)(e"*(2p3(l + w) + 93(1 - w)) + 



e-^\p3il-w) + 2q3{l + w)) 



3C cosh {2ijt) + - (3 + it;)(l + 3w) (jsge^'^* + 936"^'^*) 



(73) 
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Baryon Number Conservation Equation 



53 W 



11 9(1 -u;) 



2 2{l + w) 
(1 - w) (e'"*(2p3(l + w) + 93(1 ~ w))+ 



^^cosh i2ujt) + l{3 + w){l + iw) (pae^"* + gae-^"*) 



(74) 



where pa and qa are some constants. 



HIGHER ORDER INVARIANTS, BARYON 
NUMBER CONSERVATION EQUATION, 
GENERAL STRUCTURE 



The perfect fluid isotropic condition Eq. (fT7j) and the 
associated equation of state Eq. ([TSl) form the basic equa- 
tions for the metric components. The perfect fluid equa- 
tions of motion, on the other hand, are governed by the 
inherent energy-momentum conservation Eq. (|19l) and the 
four- velocity vector normalization condition Eq.(l5]). In 
fact, according to our discussion above those lead to the 
Euler and Baryon number conservation equations, given 
by Eg. d^ . Eg. d^ . respectively. Consequently, it is pos- 
sible to perform a consistency check, by plugging our 
solution into the right hand side of the baryon number 
conservation equation 



1 



-9 



(79) 



Higher Order Invariants 

Is the horizon surface at x — singular? To answer this 
question one may calculate some curvature scalars and 
check their regularity on the horizon. In fact, utilizing 
the MathTensor package we can work out the following 
expressions for the various curvature scalars TZ, TZ^^TZ^^, 
'R^"'^'"R^,yx„ along the metric Eq.®, 



7^ = 



24^2^(1 - iw) 

l + w 



-2,2u'^{l~iw)W(t)x + 0{x^) (75) 



576w'*C^(l -|-3w;2) 1536a;''C(l + Sw^) 



l + w 



W{t)x 
(76) 



keeping in mind that it should vanish along our solution. 
In our case of study, the line element and the four velocity 
vector are given by Eq.® and Ea. ((T^ . respectively, and 
Ea. (l79)) acquires the following form 



VRp^ cosh((p) ) + (VTr'^pT^ sinh(93) ) = 



(80) 

Plugging in now the harmonics up to third order, we 
arrive towards the following 

ipVi-i + v)W (e^ - 1 + p^-^) = (81) 



which holds thanks to Eq. (|65|) we obtained before. As 
expected, we verify that our solution respects the baryon 
number conservation equation and that the latter im- 
poses no further restrictions at the zeroth order. 



252tj4(4^(l + w) + (1 + wf +(^{9 + iiw + 9^2)) 

512a;4(2(l + w) + ^(9 + Qw + 9^^)) 
^ l + w 

where W{t) is given by 



W{t)x + 0{x^) 
(77) 



W[t) =pa (2(1 + i«)e"* + (1 - w)e-'^') + 
g3((l-He"* + 2(l + w)e-"*)- 



(78) 



As expected, those higher order invariants are regular 
on the local Killing horizon a; = and its close vicin- 
ity, signaling that there is no physical singularity in our 
solution. 



General Structure 

Let us assume that we have calculated the terms in 
our expansion Eq. (jl2p of order n — 1 and lower. Writing 
next the equation of state Eq. ([TTl) for /„ (t) and <?„ (t) and 
analyzing it, one can see it admits the following proper- 
ties: (i) doesn't comprise time derivatives of fn{t) and 
(ii) linear with respect to fn{t). Similarly to the ex- 
pression Ea.([58| this allows to represent /„(i) through 
gn{i) and its derivatives. Finally substituting the latter 
into the other Eq. (jl8l) one ends up with a linear equa- 
tion for gn{t). Carrying out this steps one may obtain 
the third order terms fsit) and gsit) in the expansion 
Eq.dlll). Those are given exphcitly by Eq.dTS]), Ea.d74l). 
respectively. 

At each order n of a; two linear differential equations 
for fiiit) and gnit) are introduced. This suggests the 
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(82) 



following double series structure 

^ oo 

T{t,x) = -x'Y.Ut)x- 

oo 

R{t,x) = 16^g„(t)a;" 

such that 

n 

rn^O 
n 

m=0 

The above structure can be put into the alternative form 

oo 
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(83) 



m— 



(84) 



R{t,x) = E bm{x)e 



m=0 



Note that the explicit Euclidean t periodicity has not 
really gone away; apart from the transformation itself, it 
enters now to the circumferential radius. 



SUMMARY 

In this paper, we have proven the existence of gen- 
eral relativistic solutions describing black holes embed- 
ded in accreting perfect fluid, and demonstrated the phe- 
nomenon for the class of P = wp equations of state. Ob- 
viously, by simply failing to respect the assumed equa- 
tion of state, the well known McVittie solution does not 
fall into this category. A perturbative analysis is per- 
formed, and explicit solutions for the corresponding met- 
ric components have been obtained. While counter intu- 
itively admitting a local time-like Killing vector on the 
event horizon itself, the various black hole configurations 
are necessarily time dependent (thereby avoiding a well 
known no-go theorem) away from the horizon. Consis- 
tently, Hawking's imaginary time periodicity is globally 
manifest on the spacetime manifold. 
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